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We calculate the phase diagram of a ladder system, with a Hubbard interaction and an interchain 
coupling t±. We use a Renormalization Group method, in a one loop expansion, introducing an 
original method to include fcii dependence of couplings. We also classify the order parameters 
corresponding to ladder instabilities. We obtain different results, depending on whether we include 
fey dependence or not. When we do so, we observe a region with large antiferromagnetic fluctuations, 
in the vicinity of small t±, followed by a superconducting region with a simultaneous divergence of 
the Spin Density Waves channel. We also investigate the effect of a non local backward interchain 
scattering : we observe, on one hand, the suppression of singlet superconductivity and of Spin 
Density Waves, and, on the other hand, the increase of Charge Density Waves and, for some values 
of t±, of triplet superconductivity. Our results eventually show that ku is an influential variable in 
the Renormalization Group flow, for this kind of systems. 
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1. INTRODUCTION 

The physics of ladder systems remains a source of con- 
siderable interest. In the last decades, many conduc- 
tors were found, with anisotropic two-leg electronic struc- 
ture, such as SVC^OsQ or Sru^ x Ca x Cu2i04 : i^ com- 
pounds. The structure of La2Cu20z\]& can also be an- 
alyzed as weakly coupled ladders, and is therefore very 
similar. The phase diagram of these compounds is very 
rich; it is now well established that these systems behave 
like Luttinger liquids at high temperature, while they can 
behave like Fermi liquids when T decreases; they exhibit 
superconducting (SC) phases of type II, which can also be 
mixed with antiferromagnetic fluctuations. In some case, 
the SC phase is found to be spin-gapped, while spinless 
phases are also reported 3 . 

From a theoretical point of view, the ladder (or 
two-coupled chain) model is the simplest quasi-one- 
dimensional one. Although all its properties are not 
entirely elucidated, it has been used by many authors 
as a toy-model, to build and explore new approaches 
(two-leaf dispersion models calculated by a Kadanoff- 
Wilson renormalization method 5| or by a bosonization 
method[(|; transition between commensurate and in- 
commensurate filling^, in particular close to the half- 
filing case0,|I| ; dimensional transition^] , etc.). These 
systems have been intensively used to inve stig ate non 
conventional SC process (with singlet 0, 0, or 
triplet 1 1 pairing). 

This paper is devoted to the study of a ladder model, 
with a Hubbard interaction (U is the Hubbard constant) 
and an interband coupling (t± is the interaction factor), 
at zero temperature. We investigate the phase diagram 
in the range of parameter < U < 87rf||?7 and suppose 
that the electronic filling r\ is incommensurate. Although 
there is yet no evidence for the existence of compounds 
lying in this range of parameter, we have hopes that some 
of these will indeed be found and confirm the theoretical 
predictions that we present here. 



The phase diagram of this system has been partially 
studied by several authors. In particular, M. Fabrizio 
have used a Renormalization Group (RG) method, in a 
two-loop expansion from the Fermi liquid solution|l2j. 
He includes interband backward scattering gb, and, 
within the range of parameter that we have investigated, 
finds a SC phase, named 'phase F, which he clearly points 
out not to be singlet s, although he did not elucidate 
the symmetry of the condensate. In his calculations, the 
RG flow of susceptibilities shows several divergences : 
the Spin Density Wave (SDW) channel coexists with the 
superconducting one. This author did alternatively cut 
the flow before the divergences and bozonize the effec- 
tive Hamiltonian; then, SDW modes disappear and are 
replaced by Charge Density Wave (CDW) instabilities. 

These results are more or less compatible with that of 
other methods, using a one-loop expansion [laJlM. °r bo- 
zonizing the Hamiltonian of the bare svstem[r7l fl8L [jjj, 
l2(il ] (see also Ref. 0), or using Quantum Monte Carlo 
method plf. These authors generally find a singlet su- 
perconducting condensate of symmetry d, which coexists 
with SDW or CDW instabilities. This complicated phase 
has also been related to the RVB phase (23] . One of the 
central question is whether the SC modes are spin-gapped 
or not, and receives various and even contradictory an- 
swers. Using a two-loop expansion, KishineJ^ observes 
a spin gap, which is suppressed when t± — > 0. This is 
also the result found with a Density Matrix RG method 
by Parkj^. 

We give a new insight into these questions, using a RG 
method, in a perturbative expansion in U. This kind of 
method has been used in many recent and very complete 
works[25l|2i|27|. Here, we calculate RG equations with 
one loop diagrams, including gb couplings, as in Ref. [l2| . 
as well as all parallel momentum dependence. 

Let us emphasize that, although we begin from the 
Fermi liquid solution, we find a phase with only SDW 
fluctuations, for small enough interband interaction tj_ < 
t± c (U), contrary to all previous results obtained by RG 
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methods, which always indicate SC instabilities as soon 
as t± ^ (see for instance Lin et al.[la^). However, this 
phase is different from the one-dimension limit. This 
is a very remarkable result, since ku is known to be an 
irrelevant variable of the RG flow 28j. However, we will 
prove in this paper that it is indeed influential in the very 
case of a ladder. 

When tj_ is increased, we observe a transition at t± c to 
a superconducting phase, where singlet SC instabilities of 
symmetry d coexist with SDW ones. This phase is found 

by many authors (see Refs. Guam and other articles 

already quoted). 

Recently, Bourbonnais et a/. |29L l3(i| have examined the 
effect of interchain Coulomb interactions for an infinite 
number of coupled chains, using RG method. Interchain 
backward scattering was found to enhance CDW fluctu- 
ations and favor triplet instead of singlet SC. We here 
investigated the effect of a Coulomb interchain backward 
interaction Cb ac k for the ladder problem. We find that 
this interaction favors triplet SC instabilities instead of 
singlet ones and CDW instabilities instead of SDW ones 
in a ladder system too. Indeed, both the singlet SC and 
SDW instabilities (if any) are suppressed when Cback is 
increased, and we observe triplet SC as well as CDW 
instabilities. The triplet SC existence region is however 
very narrow and lies inside the region t±_ > tj_ c (U). For 
large values of Cback (Cback ~ U), CDW is always domi- 
nant; this is consistent with what Lin et al. findp^f. 

On the contrary, when a Coulombian interchain for- 
ward interaction Cf or is added, all SC instabilities are 
depressed, and we only observe SDW and CDW fluctua- 
tions. 

We also present a detailed classification of the pair op- 



erators in a ladder, which are connected to the order 
parameters. It proves a very powerful tool in these so- 
phisticated RG methods. 

So, we will first give a short description of the model 
(in section 2), then present the classification of the pair 
operators (in section 3, the symmetries are explained in 
appendix 1X^1 . then we explain the RG formulation and 
techniques that are used here (in section 4). Results con- 
cerning only initially local interactions are presented in 
section 5, while those concerning the influence of addi- 
tional interchain interactions are given in section 6. In 
section 7 we conclude. 



2. MODEL 

The Hubbard model of a ladder has been studied in 
various articles. We give here a brief presentation of this 
model (see Refs. 0| or El with similar notations), 
a. Kinetic Hamiltonian 

1. The model in a 1-d representation 

The dispersion curve separates into two bands (0 and 
7r), so the Fermi surface splits into four points (— kfo, 
—kf n , kfir, kfo in the || direction, + corresponds to right 
moving particles and — to left moving ones). The bands 
are linearized around the Fermi vectors with a single 
Fermi velocity Vf (cf. Fig. We write R the right 
moving particles and L the left moving ones. Then, the 
kinetic Hamiltonian writes 



<y K K 

+ J2(K + k f0 )Ll(K)L 0a (K) + + kfr)Ll„{K)L na (Kj) . (I ) 

K K 

I 



We define the Fermi surface gap Akf = kfo — kf v . One 
then gets Akf = 2t±/vf. The discretization step in || 
direction is a, and the reciprocal vector in this direction 
is defined modulo 2ir/a. The distance between the chains 
in _L direction is b. 




Figure 1: The 2-band dispersion in || direction 



In order to give a clear representation of all instability 
processes that will be discussed afterwards, it is worth 
showing how this model can be embedded in a 2-d rep- 
resentation, which we present now. 



2. The model in a 2-d representation 

The general 2-d dispersion law writes 

e(k) = — 2£|| cos(fc||a) — 2t± cos(k±b) 

as represented in Fig. El if one writes approximately 
kf w kfo w fc/Tr, one gets Vf « 2at\\ sin(fc/a). In real 
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space, J/N corresponds to the axis along the chains, and 
y± takes only two values, y± = ±1, corresponding to 
which chain one refers to. The || axis is then discretized, 
t/ii = 1 • -N, where AN is the total number of states. 
The Fourier transformation from real space functions to 
reciprocal space ones is detailed in appendix I Dl 



From y± — ±1, one gets k± = or ±? (k± — ? and 
k± = —j are identified). Therefore, the physical states 
are limited to the four horizontal bands, shown on Fig. [21 
which are centered on each of the four Fermi points. The 
two bands centered on (±fc/ O ,0) are the left and right 
bands (L or R ), the two bands centered on (±kf 7T , j) 
are the left and right bands it (L n or R^). 



Figure 2: Two-dimensional representation of the states. The 
physical space, for a ladder, is restricted to the hatched bands 
at k± = and k± = ±?. The curves correspond to the Fermi 
surface of a 2-d system, with an infinite number of chains. The 
four Fermi points, in the ladder system, are the intersections 
of these curves with the physical bands. The symmetry C 
that maps point M onto M' is the point symmetry around 
the affine point, corresponding to fc /o+ fc /7r ^ rp^g S y mme t r y 
<P is the translation by the vector kfo — kf^, which is also 
represented by a plain arrow. 



Interaction Hamiltonian 



The most general interaction Hamiltonian can be writ- 
ten, with L k <T = LQ a (K) for k = (if, 0), L^ a = L^(K) 
for k = (K, 7r) and idem for R, 



ki.k 2 .k^,k 2 fl,tT2 
k 1 +k 2 =k' 1 +k^ 



where 5a is the two-particle coupling, and we have used 
the g-ology representation. We do not include umklapp 
interactions $-z(LLRR or RRLL) nor S4 terms (LLLL or 
RRRR) . We will alternatively use the singlet-triplet rep- 
resentation, where a takes values a = s, t, or the Charge- 
Spin representation, where a = C, S. If not necessary, we 
will omit the spin dependence a. Eqs. I@J and JJjJ give, in 
appendix I A al the usual relations between these different 
representations . 

We distinguish, following Fabrizio.l2j, g , which 
corresponds to R^L^LqRq process, g v (<-» 
RlLlLM, g f0 («-► B&LILkRq), g fn (<-► RIl\lM 



5to (<-» RlLlLirRir), g t7T (<-» B^L\L Ro), g bQ (<-> 
RqL'ILqR^) and («-> R\l\L^Rq). The definitions of 
these couplings, including the fc|| dependence are detailed 
in I A al and Fig. 1181 from symmetry considerations, one 
only needs go, g^, 17/0, gto and <7bo; m fact, we will see 
that even g^ can be deduced from go in the very case of 
a ladder, so that we only deal with four couplings. 

Bare couplings Of course, in the initial Hubbard 
model, the two-particle couplings do not depend on the 
momenta. We will define U — Ua/(jrvf) and g = 
ga/(irvf), to get rid of the physical dimensions. Then, 
, the bare couplings values are simply cji = U. 
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Additional Coulombian interchain interactions The 
above Hamiltonian H ult corresponds to local interactions. 
We have also studied the effect of additional Coulombian 
interchain interactions. 

In order to implement a backward interchain interac- 
tion, we need to modulate the bare parameters, which 
simply writes, in this case, g 01 = U + C hack , g f01 = 
U - C*back, 5toi = U- CWk and gwi = U + C hack , where 
Cback = Cbacka / ' {^v f) is the corresponding interaction 
factor. 

When we include instead a forward interchain inter- 
action, of parameter Cf or = Cf or a/(7ru/), the modu- 
lation of the bare parameters writes 502 = U + Cf or , 
.9/02 = U + Cfor, 9t02 = U - Cfor and g b02 = U - C ioY . 

Eventually, if we include both interaction, we only need 
to add the two modulations. 



c. External excitation fields 

We note Z the three-legged couplings to external ex- 
citation fields. We will write Z° w , (a = C,S) for 

Charge or Spin density waves, and Z a , (ot = s, t) 
for singlet or triplet superconducting states of symmetry 
r = s.p, d, f, g. Again, we will omit index a as soon as it 
is not necessary, and distinguish Zq W , which corresponds 
to 4i? process, z™ (<-► LlR^), z™ («-► L\R ) and 
z pw ^ LqR^). The first process corresponds to an in- 
traband mapping that relates 0-band states (horizontally 
in the 2-d representation of Fig. [2J, idem for the second 
one with 7r-band ones, while the last ones are interband 
mappings (biased in Fig. the same applies to Z sc , 
except that the processes now write LqRq, L^R^, LqR w 
and L^Rq. The definitions of all these couplings are de- 
tailed in I A bl and Fig.EO from symmetry considerations, 
one only needs zq , z^ and z+ ; again, in the very case of a 
ladder, we will see that z^ can also be deduced from zq, 
so we only deal with two couplings per instability. 

With our specific choice, the initial values of the cou- 
plings to external fields are all z% = 1. 

After this brief presentation of the model, we will now 
expound the classification of the different instabilities, 
according to their symmetries, which we will study after- 
wards. 



3. RESPONSE FUNCTIONS 

To each external excitation field corresponds a suscep- 
tibility, which is the response function of a pair operator. 
The corresponding order parameter is the mean value of 
this operator. In order to classify the different instabili- 
ties, one just needs to classify the pair operators. Their 
symmetries are detailed in appendix l"Cl 

We will first begin with SC instabilities. 



a. Superconducting instabilities 

1. SC Hamiltonian 

Let us define the superconducting order parameters 

AoV(X) = (O^(X)), where the electron -electron pair 

( r) 

operator writes O y a = 2Jxw V'x.o-IXX, X'JV'X'.ct't"^, 
with t s = io~ y for singlet states, r* x = ia x a y = —o- Zl 
r ty = io~yO-y = il and r tz — io~ z o~y — o~ x for triplet states 
{cti are the Pauli matrix, / is the 2x2 identity matrix 
and i is the imaginary number). ipjc,<r is a real space 
wave function; since electrons occupy discrete positions 
X = (ia,bj/2) (i = 1 • -N, j = ±1), we will rather write 
ipijv The corresponding Hamiltonian writes, in recipro- 
cal space variables, 

h sc = £ sc (— -<^ (t %* PlCT1 * P2ff2 

Pl:P2 

+ ^(p.P.QjfW^*^ (3) 

where Q = (Q||, Qj_) = Pi + P2 is the interaction mo- 
mentum. 

In order to simplify our expressions in this section, we 
write L p ,0,cr = L g , a (jp - k/e) and R p fi,a = Le,a(p + fc/e), 
where 9 = 0,ir. Notation f L stands for the half band 

integration J k , °_ » or IkJ 7 '- z. > which case we won't need 
to precise, since what takes place at the band limit is not 
physically relevant, in this system. 

To each order parameter A(X) corresponds an infi- 
nite number of Fourier componants, depending on the 
reciprocal space variable Q. We will only keep compo- 
nants Q = (0,0) and Q = (±Ak f , |), so that e(k) and 
e(Q — k) both lie in the physical band, close to the Fermi 
points. We will write, in short, O(0) and 0(ir±) the cor- 
responding pair operators. 

Let us now classify the different operators, according 
to their symmetry, by choosing the adequate T. The 
principles of the calculation and some details are given 
in appendix fDl 

2. Singlet s condensates 

The local pairing T(X, X') = 6^6 jj/ gives singlet con- 
densates of s symmetry. The pair operator reduces to 

O(0) component corresponds to an intraband pair- 
ing, named 0-condcnsatc (Q = 0, corresponding to Zq , 
see I A bp and writes 

O«(0) = 2^ mt ^, u +^,-it^,-ii 

i 

= J "^"-kp.oT-R-p.oj. + Lp^R-p^i 
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0{yr±) component corresponds to an interband pairing, 
named 7r-condensate (Q = (±Afcy, corresponding to 
z+, see l A b|l and writes 

i 

/adp 
-^■L p fi-[RAkf{l±l)-p,irl + Lp^ijRAk f (~l±l)- P: Ol 

+-RAfe / (-l±l)-j9,0T^P,7r| + ■^Afe/(l±l)-p,7rT^'p,0J. i 

it is however antisymmetric with parity (POi s \j) = 

-Oi s) (f )); this comes from the e~ iQ ^ b / 2 factor in the 
Fourier calculation, see details infTJl 

The s condensates are local in real space, see Fig. [3] 
(a). 

If r is replaced by 8i,i'^m$jj' , one gets extended s 
states (in reciprocal space variables, the componants are 
multiplied by cos(mPa) or some similar factor, see some 
examples in I D|) . However, we did not include these in 
our calculations. 



slightly differs); however, the signs change according to 
g symmetry. A real space representation of the different 
condensate of singlet symmetry is given in Fig. |3l 





f 

























(a) (b) (c) 



Figure 3: Real space representation of SC condensate of sin- 
glet s symmetry (a), d (b) and g (c) 

If the componants, in reciprocal variables, are multi- 
plied by cos(mPa) (or some similar factor), one gets ex- 
tended d condensates (this corresponds to the harmonic 
classification) . 

4- Triplet condensates 



3. Singlet d and g condensates 

There are also two singlet condensates of d and g sym- 
metry. 

With r = 6i t i>5j t —j/ (interchain pairing, with equal po- 
sitions on each chain), one gets another pair operator. 
O s (tt) component is zero for singlet condensate, while 
O(0) component corresponds to an intraband pairing (0- 
condensate) of d symmetry, and writes 

i 

f adp 

= J -^-Lpfl^R^pfli — L p ^-\R-p^i 
+R-p fl] L pfii - R 

With r = Sijiz$zi6j t -j> , one gets a more complicated 
pair operator. 0(tt) component corresponds to an inter- 
band pairing (7r-condensate) of g symmetry, and writes 

0(»> (tt ± ) = -2i^eT iAfc / M (^,iT^+i,-U-^,-iT^+i4i; 



■ e 2 



adp 



Ak, 



sm(a(p-k f oT—^-)) 

Z7T Z 



y^p—Akf ,7rf R±A kf-pfii + R±Ak f -p,0]Lp-Ak f ,Tri 
—Lpfi-[RAk f {l±X)-p,-K[ — RAkf{X±X)-p,it^L p fll) ; 

be careful that the symmetry of the O-condensate is 
d x 2_ y 2 (i.e. it changes sign with C, see the definition 
afterwards), while that of the 

7r-condensate is both d xy (i.e. it changes sign with 
p x and p y ) and d x 2_ y i\ moreover, d x 2_ y 2 is imper- 
fect on the 7r-condensate (for instance, it maps a factor 
sin(a(p — kfo + —j^)) onto sin(a(— p — kfo + — j^-)), which 



One also finds triplet instabilities. 

T = Si t i'^i5jf corresponds to the p x symme- 
try; o[ Pj: \o) corresponds to an intraband pairing (0- 
condensate), and writes 

O { t P *\0) = ^(^.laV'i+l.la'+^-la^+l-laOO 

= — i J {LpfloR-pfia' 

+L P -Ak f ,TTcrRAk f -p,Trc7 l ) sin(a(p — fc/o))"^' 

0[ P:c \t:±) corresponds to an interband pairing (n- 
condensate), and writes 

^ iAk,ia {^i,X^i+X,Xa' -^i,-Xa^i+X,-Xa')r^, 



+ i J 

e 2 



adp 



{Lp,OcrRAkf(X±X)-p,' 



Akf 

+ Lp-Ak f ,TrcrR±Akf-p,Ocr') Sin(a(p - kf T —^-)) T aa' I 

be careful that, because of the factor e~ l Q^i b / 2 in the 
Fourier transform, this condensate is invariant under P. 
With T = Sis> T iSj^j> , one gets an intraband pairing 

(O-condensate) of symmetry f x , given by the O t u * ; (0) 
component 

O ( /*\0) = ^ (lpi,X<r^i+X,-X<T' + A-X<T^i+X,Xa')r" a , 



i J -^j- ^ (L p -Akf,ncrRAkf- 



p,7T<T 



-Lp. 0c7 R-p, 0a ,) sixi(a(p - k f0 ))T™ a , ; 
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note that d x 2_ y 2 is again imperfect. 

With r = 5i t i>dj t -j', one gets an interband pairing (n- 

condensate) of symmetry f y , given by the Ol Jv >(iT±) com- 
ponent 

^ (TCP' 

note that p y antisymmetry is an internal one and does not 
account on the exponential factor, in the Fourier trans- 
form. A real space representation of the different con- 
densate of triplet symmetry is given in Fig. 0] 



The construction of the response functions for these in- 
stabilities is very similar to that of the superconduct- 
ing instabilities. So, we will only give the Fourier com- 
ponants of the electron-hole pair operator for Q = 
(-2fc /0 ,0), Q = (-2k f „,0) and Q = (-fc /0 - k f *, £). 

SDW and CDW operators only differ by the spin factor 
(matrix a or I), so we will also omit this factor. 



2. DW response function 

The intraband response functions write then 
r\ f ac ^P t t 

O s itc(-2fc/0, 0) = -^-R pfia L p Qa , + R p+ Ak f ,Tva L p-Ak f ,ncr' 

Obo„d(-2fc /0 ,0) = i / ^sinMe- ife ^ Q 





4 






t 








and 




f 




t 




t 




















O site (-2fc /7r ,0) 


(a) 






(b) 






(c) 







(RpfiaLp 0a ' + Rp+Akfjira^p— Akf,ncr'J 



-ikf^a 



Figure 4: Real space representation of SC condensate of 
triplet symmetry p x (a), f x (b) and f y (c) 



Extended states of the same symmetries can be ob- 
tained exactly the same way as for singlet superconduct- 
ing operators. 



b. Density wave instabilities 

1. DW Hamiltonian 

We have also investigated site density waves insta- 
bilities, defined by the order parameter A^(X) = 

Y,„AA,J>x,M*, with T ° = 1 fOT CDW > and r 8 ' = 
a x , r Sy = o y and t Sz = a z , for SDW, as well as 
bond density waves instabilities, defined by the order 
parameter Ag^X) = E^^^x+i^') 7 "™'' where 
1|| = (1,0). These couplings are intrachain, we distin- 
guish intraband and interband ones. We also investigated 
interchain couplings, defined by the order parameters 
A S T(X) = £^<<>x+g >ct <>0 ^d ADW(X) = 

E™'<<>x+g><>0> where G or G' € {1 ± ,1||+U} 
and 1_l = (0,1). 

The corresponding Hamiltonian writes, in reciprocal 
space variables, 



DW 



J2 ZP^^^UUl-Vr 



Pl,P2 
CT 1 > CT 2 



(4) 



(-^■p-Akffia^p+Akffia' + ^p^cr^p^a' ) • 

The interband response function writes 

O site {-kf - k fw , t) = i / %T-R\$o L p^a' + R t,™ L p,0a> 

J L 

_ , , . 7r f adp a k fo+ k f* 

O hond {-kf - fc^, -) = -J — e" 2 

Bin(o(p- ^Y~)){ R l ,n* L p,0a> + K+Ak f ,0a L p+Ak f • 

The above response function are intrachain ones. The 
way we have written them one just needs to add a minus 
sign before the first (or second) term of all intrachain 
operators, to obtain all interchain ones. 



4. RENORMALIZATION GROUP EQUATIONS 

a. Choice of the RG scheme 

We have used the One Particle Irreducible (OPI) 
scheme (Ref. [2^, 1^3, 0), to calculate all diagrams, 
in a one-loop expansion. We use the flow parameter 
A = A e~ e . We do not renormalize v/ nor kfo or kf v . 

b. feii dependent equations 

1. fey dependence of the couplings 

In this system, the interband backward scattering gj, 
plays a particular role. Due to momentum conservation, 
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it is not possible to put all its arguments onto the Fermi 
points. This indicates that gb is not a low energy process. 
However, it intervenes in the renormalization of low en- 
ergy processes. For instance, in the renormalization of 
go, -fr gives a contribution containing a gbo scattering, 
with a factor x+ujAfcJ' This contribution is exponen- 
tially suppressed, as A C VfAkf — 2t±. It can thus 
be neglected as soon as VfAkf is of the order or bigger 
than the initial bandwidth 2Ao On the other hand, as 
shown by Fabrizio[l2j. the gb process has to be taken into 
account if VfAkf is much smaller than A . 

Thereupon, in order to calculate the renormalization 
of gbo properly, couplings go, gf or g t with specific k\\ 

dependence are needed. For instance, ^f - gives a Peierls 
diagram : 



fc /0 -fc /0 + 2Afe f 




including coupling go(k fo ,-k fo +2Ak f ,-k fo +2Ak f ,k fo ), with 
arguments that remain separated from the Fermi points, 
even in the limit A — * 0, let us write it goi. This coupling 
separates from coupling go( fc /o,-fc/o,-fc/o,fc/o), with all ar- 
guments at the Fermi points, which we will write goo! 
therefore fen dependence is influential. This can be proved 
by comparing their renormalization, in the Cooper chan- 
nel. For instance, ^jf- gives, in the Cooper channel, a 
term proportional to g^+g^ with a constant factor, which 
is present all the way down to A — > 0. On the contrary, 
gives, in the Cooper channel, a term with a fac- 
-; the renormalization of goi 



tor 



2A 



2A+t>/|Pi+P2| 



A+vjAkf ' 

in the Cooper channel is thus exponentially suppressed, 
when A — + (for goo, the total incoming momentum is 
Pi + Pi = 0). 

We have proved that different g couplings separate 
during the flow, so the k\t dependence is hence capable 
to have an effect during the flow, when it is taken into 
account. This generalizes for g/,gt and even gb couplings. 

All this differs completely from the one-dimensional 
case, where the renormalization of the coupling with all 
momenta on the Fermi points is only governed by pro- 
cesses with momenta ±kf ± which always fall on the 
Fermi points when A — > 0. In our case, it is not possible 
to apply the same argument as soon as t± < A„. Indeed, 
we will see, in the following, that one gets different re- 
sults, depending on whether we take the k\\ dependence 
of the couplings into account or not. 

The fcii dependence can be observed, when A is large 
(and till A > VfAkf), by the separation of the different 
scattering couplings go, gfo and gto- On the contrary, if 
one puts gb = at £ = 0, this dependence is suppressed, 
and the system becomes purely one-dimensional for small 
values of t± (in that case, gb remains for all A and the 



RG equations simplify drastically, although they differ 
from the one-dimension case). 



2. fen representation of the couplings g 

In order to write explicit fen -dependent RG equations, 
one needs to define a consistent and detailed fcii represen- 
tation of the couplings g. 

Let us first note that 3 (Pi, P2, P{) corresponds 
to Pj(Pi)L^(P 2 )L(P£)R(P{), where P l are the abso- 
lute momenta in the || direction. We then define the 
relative momenta pi = Pi — kfg 1 , p 2 = P% + kfg 2 , 
P3 = P3 + kfg 3 and pi = P4 — kfg 4 , and write, corre- 
spondingly, g(pi,P2,P2iPi)- We also introduce variables 
c = Pi + P2 = P1+P2 + d, I = pi - p[ = p' 2 - p 2 + d 



and p = pi — p' 2 = p'i — p 2 + d, where d = —2Akf for 
gbo, d — 2Akf for g^, and d is zero otherwise, and then 
write, correspondingly, g(c,l,p) (d is implicit). 

At the beginning of this section, we have 
found in a diagram a particular coupling 
g (±,2Ak f ,0,2Ak f + A). when A 0, 
we get 3o(0, 2Akf, 0, 2Akt) (which also writes 
<7o(2Afcr, — 2Akf, 0) in (c,l,p) notation). Note that 
some arguments are shifted by ±2Afc/, compared to the 
coupling go(0, 0,0,0) with all momenta on the Fermi 
points. 

This could be easily generalized for all couplings g. 
So, in order to get a reasonable number of couplings, 
we have done the following approximation : all terms 
i^-, in all diagrams, are replaced by their 2Afc/ part 



(i.e. by 2Afc/ 
jer < x). 



where [x] is the biggest inte- 



_2AkfVf 

Then, it follows that we only get couplings 
<?(Pijf , 2,f , 2'P'i) ( or s( c j^p))i where all variables Pi (or 
c,l,p) are multiples of 2Afcy. 



3. fen representation of the couplings z 

All the preceding procedure generalizes to the cou- 
plings z as well. 

We first introduce a (k,c,p) representation, similarly, 
with c = pi + P2, p — P2 — Pi and k = p±, where pi 
and P2 are defined on Fig. 1191 and write, correspondingly, 
z sc (c,k) or z uw (p,k). 

Then, we apply the same approximation in order to 
get couplings, where all variables (k,c,p) are multiples 



of 2Afc 



The same conclusion applies to these couplings, prov- 
ing that their feii dependence is also influential. 



4- RG equations 

Finally, we calculate the RG flow of the 
separated following couplings : <?o(0, 0,0,0), 
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go(-2Afc/,Q,0,-2A*/), go(Q,2Akf,0,2Akf), etc. 
as well as z(0,0), z(2Akf,0), etc. 

The exact RG equations, including all || components, 
are given in I B al for the g(c,l,p), and in I B bl for the 
z sc (c,k) and z DW (p,fc). 

In order not to solve an infinite number of equations, 
we have reduced the effective bandwidth of the renormal- 
ized couplings to 4n max Afc/, where n max is an integer, by 
projecting all momenta lying out of the permitted band, 
back into it, according to a specific truncation procedure 
that will be explained after. 

We have performed our calculations with n max = 2, 3 
or 4, and the results rapidly converge as n max is increased. 



We then apply them to every coupling go, g v , etc., 
and find, altogether, exactly two independent relations 
for each one, which write, in (c, l,p) notation, 

gi(c,l,p) = gi(-c,l,p) i = 0,7r,i0,i7r ~\ 
9fn(cJ,p) = 9fo(-c,l,p) > (G-I) 

gbTr{c,l,p) = gbo(-c,l-2Akf,p-2Akf) J 



9i(c,l,p) = 9i(c,-l,p) i = 0, 7T, /0, fir ~) 
g t „(c,l,p) = gto(c,-l,p) > (G-2) 

gbnicj^p) = g b0 (c-2Akf,-l,p-2Akf) J 



5. Susceptibility equations 

To each instability corresponds a susceptibility. We 
will write Xq C ^ the different SC ones and x« W the dif- 
ferent SDW or CDW ones. 

The susceptibilities have no fcii dependence. However, 
couplings z with different ku variables appear in their RG 
equations, which we give in I B cl 

Referring to the transverse component of the interac- 
tion vectors, we will write x(0) the instabilities corre- 
sponding to intraband processes, and x(f ) those corre- 
sponding to interband ones. 

We use several symmetries, to reduce the number of 
couplings. Because of the fc|| dependence, it is not as 
easy to apply them as in ordinary cases. We give here 
some indications, which are completed in appendix. 

c. Symmetries 

1. Ordinary symmetries 

We apply C the conjugation symmetry (C : r — > 
r, p — > p, a — > —a), A the (antisymmetrical) exchange 
between incoming particles, A' the exchange between 
outgoing particles, P the space parity (P : r — > r, p — » 
— p, a — ► a) and S the spin rotation (S : a — > —a). We 
will also use p x , p y (the mirror symmetries in the || and 
_L directions), f x and f y . Note that P — p x p y - 

H cin (Eq. (jj) and H int (Eq. J5J|) satisfy all these 
symmetries, whereas SC instabilities, governed by -ffsc 
(Eq. (5)). are not invariant under S or C, which allows a 
natural classification of the states, and DW instabilities, 
governed by H BW (Eq. @), do satisfy CS, AS or A'S, 
but not C, A, A' nor S" symmetries. 

All the relations satisfied by S or Z couplings are de- 
tailed in I C al From what precedes, one will not be sur- 
prised that those for 9 couplings are simpler and less 
sophisticated than those for Z ones. 

Relations of couplings g We are not interested here 
in the symmetries that relate, for instance, a LRLR cou- 
pling to a RLLR one. Instead, we only keep RLLR cou- 
plings and deduce all the symmetries that keep this order. 



One observes that G-f relates g/o to gf n and g^ to 
gbn, while G-2 relates gto to g^ and gbo to gbn- The 
combination of G-1 and G-2 thus relates gfo to gf n and 
gto to g t7T - 

Relations of couplings z We similarly deduce all sym- 
metries that keep the LR order; we apply them to every 
coupling zq, Zn, etc., and find only one relation for each 
one, which writes, in (fc,c,p) notation, 

ii {c,k) = z g(- c ' A; - c ) * = °>*] 

«f°(c,*) = ±zf°(-c,jfc-c) f ^- bL - 1 > 

zfS(c,k) - ± 2t f (-ck-c) J 

where ± reads + for r = s or T — p x and — for Y = g or 
T = f y . Note that 2q or s^ c correspond to intraband 
condensates, while or s 530 correspond to interband 
ones. 

zf w (p, k) = ±z 9 DW (p, -p-c) e = 0,n\ 
z DW {p,k) = ±z^ w {p,-p-c) J 

(Z-DW-1) 

where ± reads + for site ordering and — for bond order- 
ing. 

One observes that Z-l relates z + to z_. 



2. Supplementary symmetry 

As we already noted, ordinary symmetries do not relate 
go to g v , nor Zq to z v . However, since we choose identical 
bare values at I — 0, and since the RG equations are 
symmetrical, we observe an effective symmetry between 
these couplings : we will show here that this does not 
occur by chance, but that it follows a specific symmetry 
C, which only applies to the ladder system. 

C is a kind of conjugation : it generalizes the electron- 
hole symmetry as follows. 

Let us first consider the case of a single band one- 
dimensional system; we find an electron-hole symmetry, 
described in Fig. (a). 
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k f +p 



(a) 




Figure 5: Symmetry around the Fermi points : 
system; (b) in a 2-band system 



(a) in a 1-band 



where ± reads + for r = s or T — p x and — for r — d, g 
or T = f x ,f y . 



z™(p,k) = ±z™(-p,-k) 
z^ w (p,k) = ±z™(-p,-k) 



(Z-DW-2) 



where ± reads + for site ordering and — for bond order- 
ing. 

Again, one verifies that zq and related (as well 

as z + and zJ). 



For R particles, it conjugates an electron with momen- 
tum kf + p and a hole with momentum kf — p. In the 
momentum space, it is a symmetry around kf. One can 
write Cipp — V4fe -p and Ctpj, — ip2k f -p- For L particles, 
the same relation applies if one simply changes kf into 
-k f . 

Let us now return to the two-band system. This sym- 
metry generalizes by turning the momenta around the 
isobarycenter of the Fermi points, as shown in || space in 
Fig. 03(b). For R particles, (kf + fcj 7r )/2 points to the 
isobarycenter and you now get Cip p — ^kto+kt -p e ^ c - 

In the two-dimensional representation, C-symmetry is 
a point symmetry around (± f0 ~^ f " , jr), as shown on 
Fig. [3 (the sign depends on whether it is a R or L mo- 
mentum) . 

Actually, there is an alternative symmetry, which we 
write fD, which also maps the band k± = onto the band 
ir : it is a translation by the vector (±(fc/ 7r — fc/o), ?), as 
shown on Fig. [3 Some of the bare couplings satisfy this 
ft symmetry (g , gf,g t , z or z*), but some don't (g b , z + , 
Z-). Since interactions are mixing all couplings as soon 
as the flow parameter £ > 0, the renormalized couplings 
will break the ft symmetry. Therefore, we can't use it. 

On the contrary, one verifies that H c [ n , H lrAl Hsc and 
flow are invariant under C. The induced relations sat- 
isfied by S or Z couplings are detailed in I C bl 

In fact, C and ft weakly correspond to the d x 2_ y 2 sym- 
metry in two dimensions. 

Supplementary relation of couplings g It is straight- 
forward that C keeps the RLLR order when one applies 
it to any coupling gi\ so we find a new relation for each 
one, which writes, in (c, l,p) notation, 



Truncation 



gfn(c,l,p) 
9tm(c,l,p) 



g (-c,-l,-p) ^ 
9fo(-c,-l,-p) 
9to(-c,-l,-p) 
9w(-c, -I, -p) . 



(G-3) 



One verifies that go and g w are related; in fact, G-3 relates 
all O-couplings to 7r-couplings. 

Supplementary relation of couplings z Similarly, C 
keeps LR order when we apply it to any coupling Zi] 
so we find a new relation for each one, which writes, in 
(k,c,p) notations, 



r sc 



c, k) 



z s _ c (c,k) = ±z s + c ( 



±4° 



-k) 
-k) 



Understanding symmetry relations does not only help 
us to reduce drastically the number of couplings, it is also 
an essential tool to make a proper truncation procedure, 
as we will explain now. 



Triplet notation 



de- 



(Z-SC-2) 



Let us first introduce a useful notation for the Ay 
pendence. 

We have already defined the relative momentum repre- 
sentation g{pi,p%, p'2, Pi), as well as the g(c, l,p) notation, 
and explained how to keep only couplings, the arguments 
of which are all multiples of 2Afc/. We will focus on the 
(c,l,p) notation and write c = 2niAfc/, I = 2ri2Afcy, 
p = 2n^Akf, with (711,712,713) G Z 3 . 

In short, we can write <?i(ni, n2, v.3) (i — 0, /0, £0,60), 
where (711,712,713) is called a triplet. Mind that, us- 
ing symmetry relations, two triplets (711,772,713) and 
(71^,712,713) can represent the same coupling. One says 
that these triplets belong to the same symmetry orbit 
(or symmetry class). Mind also that the orbits are dif- 
ferent for each coupling go, gfo, gto and gbo- 

It would take too long to give an exhausted list of these 
orbits. Let us just observe that (0,0, 0)'s orbit has only 
one element (itself), except for 175, the orbits of which we 
detailed in appendix I C cl 



2. Choice of the truncation procedure 

Obviously, one needs only renormalize one coupling 
per orbit. From the fundamental rules, explained in 
14 b 21 it follows that, even if one starts with only 
^(0,0,0), 3/0(0, 0,0) and 5*0 (0,0,0), the RG equations 
will generate an infinite number of orbits. So, we will 
only keep couplings which satisfy \rii\ < ro max (for a given 
Timax); but even so, in the RG equations of some orbits 
intervene couplings, with arguments lying outside of the 
permitted band (i.e. the distance of the corresponding 
momentum to the Fermi point exceeds 2n max Afc/). In 
order to get a consistent closed set of differential equa- 
tions, one needs to put these extra couplings back, inside 
the set of allowed couplings. 



For instance, imagine that n max = 2, and that g(3, 2, 2) 
coupling intervenes in a RG equation. One cannot, un- 
fortunately, just map (3,2,2) i— > (2,2,2), because these 
triplets do not belong to the same symmetry orbit. In do- 
ing so, one would get a very poor k\\ dependence; we have 
actually proved, in the case of a (pi,P2,P2) representa- 
tion, that all orbits of go except (0, 0, 0) would collapse 
into one single orbit. 

Therefore, the truncation procedure must be compati- 
ble with all the symmetries of the system. We have used 
the (c,l,p) notation, which is very convenient. One can 
check that all the symmetries conserve c + I + p modulo 
4Akf. This explains why ni i— > n, ± 1 can't be compati- 
ble with the symmetries. On the contrary, rii i— > n., ± 2 is 
completely compatible, i.e. it maps a triplet onto on an 
already defined orbit; hence we have used this mapping 
for the extra couplings. 
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Figure 6: Flow of the susceptibilities xf ^ an d Xs W (intra 
or interband), at U = 1 and 2t±/Ao = 1.94 : (a) usual RG 
procedure; (b) including k\\ dependence. You observe that 
Xs (0) does not diverge but only reaches a plateau. 



With n max = 2, for each coupling g we find 63 different 
couplings 1 , which separate into 23 orbits (having 1, 2 or 
4 elements), except for the gb couplings. For these, the 
enumeration is more tedious, we eventually find 8 orbits 
(of 4 elements, see l C c|) . There are altogether 3 x 23 + 8 = 
77 different orbits; if we include the spin separation, we 
thus need to calculate 154 coupled differential equations. 
n max = 3 gives 390, while n max = 4 gives 806. 



e. Divergences of the susceptibilities 



In the range of values for U that we have investigated, 
the RG flow is always diverging. 

When the initial interaction Hamiltonian iJi nt is purely 
local, i.e. when the interchain scattering are discarded 
(Cback = Cf or — 0), the interband SDW susceptibility 
Xs W (f') ' s always divergent. In the superconducting 
phase, the SC singlet d susceptibility xf C ^(0) 1S also 
divergent, at the same critical scale A c = Aoe^ c . A 
third susceptibility Xs W (0) increases and reaches a high 
plateau (see Fig. 0) . Almost all other susceptibilities re- 
main negligible. 



When the parameters Cback or Cf or are increased, 
both x<j W (§) and xf ^ decrease; they are progressively 
replaced by the divergence of the CDW susceptibility 
Xc (f )> an d' m t ne case °f backward scattering, of the 
triplet SC susceptibility xf C(/) (0)- 

One finds at most four divergent susceptibilities 
(X s s° (d) (0), xg W (-F ), X t SC(/) (0) and xg W (f )) at a time. 

Since the RG flow is diverging, we cannot further cal- 
culate the renormalized couplings. To deduce a phase 
diagram, we must find out which mechanism dominates; 
we used two different criteria : according to the first one, 
we simply take the susceptibility which reaches the high- 
est value |x| at A c ; according to the second one, we take 
the susceptibility which has the highest slope. 

These two criteria bring non equal results. Although 
the first one is a poorer criterion, its conclusions remain 
stable when either the precision or n max are changed. 
The second is however preferred, as we will see its con- 
clusion are physically consistent, contrary to the first one. 



PHASE DIAGRAM WITH INITIALLY 
LOCAL INTERACTIONS 



Let us first discuss the case of initially local interac- 
tions (Cback = Cfor = 0, no interchain scattering). Of 
course, we can only fixe Hint at £ = 0, and the flow will 
develop non local interactions. 



Results 



1 that are {(i, j, k), i, j, k=-2, 0, 2}(J{(±1, ±1, i), i=-2, 0, 2} 
U{(±1, i, ±1), i=-2, 0, 2}{J{(i, ±1, ±1), t=-2, 0, 2} 



We begin with the phase diagram obtained when the 
k\\ dependence is neglected. 
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Figure 7: Flow of the couplings goi and <?02 at : (a) t± < t± c (U); (b) t± = t± c (U); (c) t± > t± c (U) 



Phase diagram with no kn depende 



In the region of the phase diagram that we have in- 
vestigated (0 < U < 2), the SC susceptibility x sc(d) (0) 



as well as the SDW susceptibility 



is always divergent 
Xs W (f) (see Fig.©. 



According to the slope criterion, xf ^-* (0) always dom- 
inates. We induce that this region is superconducting 
(this is consistent with the conclusions of Fabriziop"^). 
and that the pairing is of symmetry d; however, the pres- 
ence of SDW instabilities, developing in the same region, 
makes a detailed determination of this phase very uneasy 
and beyond the possibilities of our approach. We will call 
it SC phase. 




-Xf (d) (0) 
-Xf c(,) (0) 

DW 

<°) site 

DW 
xXg (0) site 

»X c DW (l)si,e 

< W (l)si,e 

^ W ©bond 

< W (g>bond 



Figure 8: Flow of some susceptibilities, at U = 1 and 
2tx/Ao = 0.016. xf C ' d ' ) (0) is almost suppressed at A c . 



2. Phase diagram with k\\ dependence 



SC phase When t± is large, we find similar results. 
For instance, with [7 = 1, there is no significant differ- 
ences for 2t ± /A > 1.94 (see Fig.©. 

SDW phase On the contrary, the superconducting 
susceptibility is almost suppressed when t± is small 
enough. For U = 1 and 2t ± /A = 0.016, x S s C(d) {0) is 
5 orders of magnitude smaller than Xs W (f )> at A. c , see 
on Fig. 121 In this phase, the SDW 's instability develops 
so rapidly that it overwhelms all other processes. This 
indicates indeed a pure SDW phase. 



Hence, this result differs drastically from those ob- 
tained when the fcii dependence is neglected. 

Moreover, we observe a transition between the SDW 
phase and the SC one, at a critical parameter t± c (U). 

Critical behavior We characterize this transition by 
different ways. 

First of all, the behaviour of the renormalized two- 
particle couplings change very rapidly, at t± c (U): as t± 
decreases, |<?oi(A c )| and |<?02(A C )| shrink suddenly, then, 
after a little interval, <7oi(A c ) and 502 (A c ) become pos- 
itive (and even > U, see on Fig. [7J). We also observe 
changes, though less significant, for the other couplings 
g (for instance, |p/oi(A c )| decreases and <7&i(A c ) becomes 
negative when t± increases). 

Moreover, we observe (on Fig. © a marked site/bond 
separation of the SDW susceptibility, at t± c (U). The 
site and bond SDW susceptibilities are degenerate for 
t_L < t± c {U) (in the SDW phase), which is consistent 
with the fact they should be equal at tj_ = (where the 
system is purely one-dimensional, see Ref. |23), while 
they smoothly separate after the transition. The same 
site/bond separation occurs for the CDW susceptibility. 
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Figure 9: Curves of xf c(d) ( A c) and x<? W (A0 versus 2i x /Ao- 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



— '±C2 


(n=2) 




1 1 1 


. — '±C 


(n=2) 






- '±C2 


(n=3) 




/ / 


t ±c 


(n=3) 






— '±C2 


(n=4) 




/ / 


... t ±c 


(n=4) 


SC 










transition 
























SDW 






1 1 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
U 

Figure 10: Phase diagram when the ku dependence is in- 
cluded. 



Transition region The behaviour of most of the pa- 
rameters that we have examined indicates the same crit- 
ical value t± c (U), which we have determined exactly, us- 
ing the slope criterion. 



However, |xf C(d) (0)1 < IXs W (f)l holds until i ± reaches 
a value t± C 2(U); this second critical value, which cor- 
responds to the height criterion, is confirmed by minor 
modifications of behaviour, which occur in the interval 
t± c (U) < t± < t± C 2(U) and are very smooth (for in- 
stance, <?&i(A c ) and #b2(A c ) cross) . 



The numerical determination of t± C 2(U) is very stable 
(see Fig ll(J[) . and the complete behaviour, from the SDW 
phase to the SC phase, is clear on Fig. which shows 
the absolute values of the susceptibilities at A c . 



The region t± c (U) < t± < t± C 2(U) is called transition 
region, we believe it is a superconducting phase, where 
SDW instabilities seem however to dominate. The SDW 
are precursor manifestations of the pure SDW phase, 
which is next to the transition region. 



As already stated, our results converge very rapidly 
when the band width on which we project the momenta 
is increased. The complete phase diagram (for Cb ac k = 
Cfor = 0) is shown on Fig. rni| for n max = 2, 3 and 4. 



Moreover, it is most interesting to note that, contrary 
to the pure SDW phase, the transition region can be de- 
tected when the A;|| dependence is neglected (the value of 
t±c2(U) is lowered). This can be seen, for instance, on 
Fig. fT^ which corresponds to Fig. [5] with no k\\ depen- 
dence. 
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Figure 11: Curves of xi d) (A0 and Xs W ( A J 

versus 2t±/Ao 

when ku is neglected. 



SC critical temperature The critical energy A c at 
which SC susceptibility diverges gives an approximate 
indication of the SC critical temperature. We give a 
plot of A c /Ao versus 2tj_/Ao; as seen on Fig. A c is 
roughly decreasing with t±. The band gap parameter 
tx is also increasing with pressure; therefore, this be- 
haviour is compatible with the experimental data, which 
show that T c is decreasing with pressure in quasi-one- 
dimensional organic compounds. 
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Figure 12: Curves of A c /Ao = e c versus 2fx/Ao with 
n max = 2 (a), n max = 3 (b) and n max = 4 (c). 



b. Discussion 

As already mentioned, as soon as the dependence of fcii 
is included, we observe two separated phases, one purely 
SDW, the other one a SC phase with competing SDW 
instabilities. 

On the other hand, for fj_ > 1, i.e. when the initial 
bandwidth lies inside [kf n , fc/o], the fc|| dependence has 
no observable influence on the susceptibilities. 

In the SDW phase, our results prove the existence 
of large antiferromagnetic fluctuations. We believe that 
these SDW instabilities are not the signature of a local- 
ized antiferromagnetic ground state, but of antiferromag- 
netic itinerant electrons, as it is indeed observed in Bech- 
gaard salts. Actually, the flow is driving towards a fixed 
point, which does not seem to be the one-dimensional so- 
lution : for instance, the renormalized couplings gi and 
<?2 of the 1-d solution are and 1/2 and differ from the 
values which we obtain when the flow is diverging, in the 
SDW phase (see Fig. □ (a) ) . 

We induce that the spin-gap should disappear in this 
SDW region, which is consistent with what Park and 
Kishine|23 l l24 | claim. 

In the SC phase, the SC divergence is due to the 
Cooper channel, while that of density waves is due to 
the Peierls one (see, in the case of a single band model, 
Refs. [13, EfJ). The appearance of d-wave supercon- 
ductivity in ladder systems is well understood within a 
strong coupling scenario, where a spin gap leads to in- 
terchain Cooper pairing. However, in our calculations, 
we see that superconducting correlations are always en- 
hanced by SDW fluctuations. Contrary to what Lee et 
al. claimed first |37j. there is an itinerant electron mech- 
anism in this case, which is the weak coupling equiva- 
lent of the localized electron mechanism in the strong 
electron scenario. It was proposed by Emery [3^, and is 
essentially the spin analog of Kohn-Luttinger supercon- 
ductivity. The mutual enhancement of the two channels 
is also discussed in Refs. |33l l39l|. 



As a consequence of this mutual enhancement, the 
spin-gap should not appear with the first appearance of 
SC instabilities, but for somehow larger values of t±. 

Moreover, we observe that the SC pairing is a Q = 
mechanism, while the SDW are excited by Q = 
(2Akf, vectors. This can be explained by the symme- 
try of each channel. The Green function of the Cooper 
channel gives a factor l/(e(k) + e(— k)) and is minimized 
with the k >— > — k symmetry, which corresponds to an in- 
traband process. The Green function of the Peierls chan- 
nel gives a factor l/(e(k) + e(k + Q)) and is minimized 
with the k i— > k + Q symmetry, which corresponds to an 
interband process. 

This can also be seen in the RG equations. 
d\n(z sc (0))/dl depends only on g and g t , whereas 
d\n(z sc (j))/d£ depends on 5/ and gb- Since gb processes 
are depressed as soon as A < 2t±, 0-condensate are fa- 
vored. Moreover, this predominance is stabilized by the 
dg t /d£ equations, in which the Cooper term depends on 
go, and by the dgbjdt equations, in which the Cooper 
term depends on j/. 

The same argument applies for DW instabili- 
ties. dln(z DW (0)) / dl depends on go and gb, whereas 
d\n(z BW (?)) / dl depends on gj and g t . So, 7r-processes 
are favored. Again, this is stabilized by the dgt/dl equa- 
tions, in which the Peierls term depends on gf, and by 
the dgb/di equations, in which the Peierls term depends 
on go- 

The critical temperature, in the SC phases, is indi- 
cated in Fig. We chose U = 0.6 in order to avoid 
the SDW phase. The general trend is that of a quasi- 
one-dimensional system; the increasing curve, for small 
values of t± , can be related to transition effect and to the 
furthered influence of the SDW fluctuations. 



6. PHASE DIAGRAM WITH INITIAL 
COULOMBIAN INTERCHAIN SCATTERING 

a. Results 

1. Influence of a backward interchain scattering 

Let us now study the effect of a backward interchain 
scattering. This type of coupling has been inverstigated 
by Bourbonnais et aZ.[2£| in the context of correlated 
quasi-onde-dimensional metals, for which CDW correla- 
tions are enhanced and triplet superconducting instabil- 
ities can occur. When the parameter Cback is increased, 
the behaviour of the susceptibilities depends on the pa- 
rameters (ij_, U). 

Appearance of triplet SC and CDW For ij_ < t± c (U), 
the SDW phase exists for Cback small enough. As Cback 
is further increased, the SDW instabilities are replaced 
by CDW ones. The transition is smooth, and there is a 
narrow region where both SDW and CDW coexist (region 
@ in Fig. EE). We show a 2t ± /A = 0.01 section of the 
susceptibilities at A c on Fig. ^| 
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For t± > t± c (U), the SC phase (with SC singlet d 
and SDW instabilities) exists for Cback small enough. As 
Cback is further increased, the singlet SC modes are re- 
placed by triplet ones, while SDW are replaced by CDW. 
Singlet and triplet SC appear to be antagonistic, and the 
transition is very pronounced; in the coexistence line be- 
tween them, one also finds SDW and CDW divergences 
(see Fig. I14|) . On the contrary, the transition between 
SDW and CDW is very smooth, although the coexis- 
tence region is still narrow (region @ in Fig. I16[) . We 
show a 2ij_/Ao = 0.1 section of the susceptibilities at A c 
on Fig. US 



When Cback is large enough, the triplet SC modes are 
suppressed, and the region is a pure CDW phase. We 
show a section of the susceptibilities at A c , for Cback = 
0.2 and Cback = 0.3, on Fig. El which clearly indicates 
the domain of existence of the triplet SC. 




( SC(d) ( 

wrSC(f), f 



'(0) 



t xg w '°' site 

*Xs>)site 
» Xg w (g) site 
°Xg W (g)site 
■ Xg W (g)bond 
o X° W (g) bond 



Site/bond separation As we have already observed it, 
in the case of Cback = Cf or = 0, for small values of fj_, site 
and bond SDW susceptibilities are degenerate, as well as 
site and bond CDW ones. 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 14: Flow of the susceptibilities for 2t±/Ao = 0.32 and 
Cback = 0.18. 

The triplet SC condensate has f x symmetry. The cor- 
responding susceptibility is mostly divergent in a region 
of coexistence with SDW and CDW (region © in Fig. lTfi|) . 
but it is also divergent in a region of coexistence with only 
CDW (region (T) in Fig. OH). 



This generalizes for all values of Cback- The site/bond 
separation line is an increasing function t± g (Cback) of 
Cback, shown on Fig. ^| For small values of U, this line 
crosses the SC domain, but for U =1 it is already discon- 
nected from the SC frontier (although it remains close to 
it). 
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U=1 



2ti/A„ . 
1.28- 




0.2 0.4 0.6 0.8 1.0 Miack 

Figure 16: Phase diagram for [7=1. The shaded area indi- 
cates the divergence of the triplet susceptibility. The dashed 
line separates site/bond degenerated states (below) and non 
degenerated ones (above). Other lines and domains are ex- 
plained in the legend or in the text. 



2. Influence of a forward interchain scattering 



The phase diagram when Cf or is included is very rich, 
and beyond the scope of this article. 

We would like to emphasize only the fact that all SC in- 
stabilities are suppressed when Cf or is increased. Fig. El 
gives a typical flow of the susceptibilities, with a large 

Cfor • 




+ x s (0) interchain 
x X™(0) intrachain 
» X° W (b' in < ercn ain 
H Xo W (l) intrachain 



° X s (5) intrachain 



Figure 17: Flow of the susceptibilities, t± = 0.02 and : Cfor 
0.95. 



b. Discussion 

Let us analyze these behaviors, which follow simple 
trends. 

The CDW instabilities are enhanced when gc is in- 
creased, whereas SDW ones are enhanced when gs is in- 
creased (this can be verified in the corresponding RG 
equations of IB b|l . Similarly, singlet SC instabilities are 
enhanced when g s is increased, whereas triplet ones are 
enhanced when g t is increased. 

So, an increase of Cb ac k implies an increase of the real 
space coupling g±, and thus, from Eq. JSJ, it favors CDW 
instabilities against SDW ones, and from Eq. Q, it favors 
triplet SC instabilities and depresses singlet SC ones. 

Similarly, an increase of Cf or implies an increase of the 
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real space coupling 52, and thus, from Eq. ijSJ, it favors 
CDW and SDW instabilities, and from Eq. Q, it de- 
presses SC ones. 

Of course, we examine here the influence of parameters 
Cback an d Cback on the bare couplings. However, we 
believe that the flow could not just simply reverse this 
influence, even if the renormalized values of the couplings 
differ a lot from their bare values. Moreover, one verifies 
that these conclusions exactly correspond to the observed 
behaviour. 

The density wave interactions are on site, whereas the 
SC pairing are inter-site (except for singlet s one), so we 
believe that the DW instabilities appear first, and then 
enhance the SC ones. This is not true of the DW bond 
correlations, but we observed no divergences of these 
ones, and we have not studied any other sophisticated 
inter-site DW excitation response. 

From this point of view, the fact that ir DW processes 
are favored, as we already discussed before, implies a 
7r dephasing between both chains of the ladder. The ir 
dephasing of the SDW thus fits perfectly singlet d con- 
densate (which consists in a pairing of two electrons on 
a rung, with opposite spins, see Fig. |3] (b)). This ac- 
counts nicely for the appearance of singlet SC instability, 
induced by SDW one. 

In the same trend of ideas, the ir dephasing of the CDW 
fits triplet f x condensate (which consists in a pairing of 
two electrons on each chain, one stepped by unity from 
the other, see Fig.0](a)) and accounts for the appearance 
of triplet SC instability, induced by CDW one. 

On the contrary, the triplet p x condensate consists in 
two following electrons on one chain (see Fig.^(b)), this 
pairing is not enhanced by CDW instabilities; in fact, it 
is the analog of singlet s condensate, which is not either 
enhanced by SDW instabilities, and is therefore disad- 
vantaged, compared to d pairing. 

As can be observed on Fig. triplet f x condensate are 
not incompatible with CDW. For instance, one could eas- 
ily figure out a succession of condensate, with alternate 
spins, inducing back a global modulation of the chains. A 
similar scenario is not possible with triplet p condensate. 

One should be aware that the symmetry classification 
we have used is very specific of the ladder system, and 
could not be extended to an infinite number of chains. 
The difference between p and / condensate is very sub- 
tle and the situation could reveal quite different in the 
general quasi-one-dimensional systems. 



7. CONCLUSION 

We have investigated the phase diagram of a lad- 
der system, in the Hubbard model, with an interchain 
coupling tj_, using functional RG method, in the OPI 
scheme. We have introduced an original parameteriza- 
tion of the fc|| dependence, and obtained rather new re- 
sults, in particular, we have proved the existence of a 
new phase with only SDW fluctuations, for small enough 



values of t± . From the divergences of the scattering cou- 
plings, we induce that this phase is different from the 
one-dimensional solution. However, for very small values 
of t± (t± < AolO~ 4 ), we find the usual one-dimensional 
behaviour. 

Our results altogether prove that the k\\ dependence 
is important and must be taken into account in such a 
ladder system. The fact that this variables become in- 
fluential in a ladder does not mean that the correspond- 
ing couplings g(2Akf, — 2Akf, 0), etc., are relevant. In 
fact, if the cut-off A — > 0, these couplings are left out of 
the integrated band, so they could only be marginal j2^. 
However, the divergence takes place at A c , which is of 
the order of Akf, and this explains why these couplings, 
which arc shifted by ±2 Akf from the Fermi points, have 
a non trivial behaviour and have to be taken into account. 
Moreover, as already explained in 14 b II during the flow, 
these couplings influence those, with all momenta at the 
Fermi points, until A = 2Afc/. This influence is still sen- 
sitive, when the divergence takes place. This explains 
why we could distinguish a new phase, which has not yet 
been observed by usual methods. 

When t±_ is very large, however (for instance, t± ~ Ao), 
the flow continues up to A c <C VfAkf (otherwise, the 
integrated band would not vary much and the renormal- 
ized couplings neither differ much from their bare values), 
and the above argument applies, proving that couplings 
g(2Akf, — 2Akf, 0), etc., are marginal or irrelevant. In 
that case, k\\ are not influential, and our results coincide 
indeed with former calculations. 

We have also given a detailed classification of the re- 
sponse function, which provides a convenient tool for the 
determination of order parameters and of related suscep- 
tibilities, corresponding to different instability processes. 

We are proceeding now to a complete study of the long 
range correlations, and in particular, of the uniform sus- 
ceptibility. This task however proves quite difficult, be- 
cause of the dependence, which has to be carefully 
taken into account. We expect that the spin-gap will in- 
deed disappear in the SDW phase we have brought to 
evidence. 

We have also investigated the influence of interchain 
scattering, and showed that a backward interchain scat- 
tering can raise triplet superconductivity, a resulst con- 
sistent with the conclusions of a previous work by Bour- 
bonnais et aZ. |29l . l30l| on correlated quasi-one-dimensional 
metals. The appearance of triplet SC in a ladder 
is a very exciting and promising result, since various 
authors[43, claim to have found experimental evi- 
dence of these instabilities. Even the narrowness of the 
triplet SC existence region seems to fit the experimental 
data, which report high sensitivity of these fluctuations 
to some key parameters. This work gains to be compared 
with the previous work of Varma et al., who did similar 
investigations [42^ . 

We would like to thank N. Dupuis, S. Haddad and B. 
Dougot for fruitful discussions and advices. J. C. Nickel 
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wishes to thank the Gottlieb Daimler- und Karl Benz- 
Stiftung for partial support. 



Appendix 
A. Couplings 

a. Two-particle couplings 

Here are the definitions of the different couplings g, 
from the two-particle parameter 9, and the corresponding 
diagrams. 







9o{Pl,P2,P' 2 ,Pl) = 




S(fc /0 


+ Pl, 


-kfo+P2, -kf +p' 2 ,kfo - 

9iv(Pl,P2,P 2 ,p'l) = 


"Pi) 




+ Pl, 


-kf„ +P2, -kfr + p 2 ,k f7r - 

9f0(Pl,P2,P 2 ,Pl) = 


fPi) 


S(fc/o 




-kfir +P2-, -kfn +p' 2 ,k f0 - 
gfir(Pl,P2,P 2 ,Pl) = 


tVl) 












-kf + P2-, -kf + p' 2 ,k f7r - 

9to(Pl,P2,P2,Pl) = 


>rP'l) 


S(fc/o 


+ Pl, 


-kfo + p 2 , -k f7r + p' 2l k f7T - 

9t7r(Pl,P 2 ,P 2 ,p'l) = 


>rPl) 




+ Pl> 


-kf w +p 2 , -k f0 +p' 2 ,k f0 - 

gbo(pi,p 2 ,p' 2 ,p'i) = 


^Pl) 


S(fc/o 


+ Pl, 


-fc/7r +p 2 , -kf +p' 2 ,k f K - 

9b7r(Pl,P 2 ,P 2 ,p' 1 ) = 






+ Pt, 


-kfo+p 2 , -kfr +p' 2 ,kf - 


T-p'l) 



.90 



9to 



.960 



X 

fc /0+ p l ~ k fir JrP 2 

X 

k fO+Pl ~ k f-7v+P2 



~ k fO+P2 fc /77-+P'l 

X 



.9/7, 



.96 jt 



X , 
X , 

fc /tt- +P1 — k fO+P2 

X , 

~ fc /7T +P2 k fO+Pl 

fc +P1 ~ k fTV +^2 

X 

-fc/0+P2 fc /0+Pl 



- k f 7v +P2 k f^+Pl 

Figure 18: Schematic definitions of the couplings g 



b. Other couplings 

Here are the definitions of the different couplings z, 
from the couplings to external fields Z. 

We omit the spin index a nor the symmetry index T, 
and use the notation explained further in appendix I C al 
Mind that T — s, d, g for a = s (singlet) and T = p, f 
for a = t (triplet). The symmetry (s, d x 2_ y 2, g, p x , f x 
or f y ) applying to each one is detailed in the main text. 



zq (pi,p 2 ,q) = 

+ Px,kfo +P2, (g,0),0) 
z^ c (pi,P2,q) = 

Z sc (-k f7r +pi,k fn +p 2 ,(<Z,0),7r) 

z+ c (pi,p 2 ,q) = 

Z sc (-k f ^ +pi,k f0 +p 2 , (q, f ),0) 

z sc (p 1 ,p 2 ,q) = 
Z sc (-fc /0 + pi, */„ + pa, (q, f ), tt) 

za W (pi,P2,q) = 
Z DW (-fc /0 + Pi, k f0 +P2, (q - 2fc /0 , 0), 0) 

z^ w (pi,p 2 ,q) = 

Z BW (-k f7T + Pl ,k f7! + p 2 , (q - 2fc /7r , 0), tt) 

z+ w (pi,P2,q) = 
Z DW (-fc /7r +pi, kf +P2, (q - kf - kfn, f ),0) 

z DW (pi,p 2 ,g) = 
Z DW (-fc /0 +pi, fc /7r +P2, (<? - fc/o - f ),tt) 



2 



•X. 

fc 

•X 



DW 

z 



9~ 2fc /0 



k fO~ k fir \ 



k fO+P2 

- k fa+pi 

k fO+P2 

fc /0+P2 
" fc /0+Pl 

" fc /7r +P1 



•X. 

fc 

k 



fc /tt +P2 

- fe /0+Pl 
fc /tt +P2 

~ fc /7T+Pl 

fc /tt +P2 
-fc/O+Pl 



Figure 19: Schematic definitions of the couplings z 



The relations between the different representations can 
be found in Refs. [2]j or 01 . Here, they reduce to : 



9s = -gi - .92 / 4 n ,9c = .92 - 2#i 
.9* = .9i ^ 92 9s = .92 



(5) 



B. RG equations 

We give here the detailed RG equations. 
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a. g couplings 




Here are the RG equations for the couplings g, in 
(c,l,p) representation. 

The spin dependence is given, for all terms, by 

0,1 

where C and T correspond, respectively, to the Cooper 
and Peierls channels, and are given, in the g-ology repre- 
sentation, by 

;) *--(; :) 

i 



see, for instance, Refs. [13 E3 E3 . In the following equa- 
tions, all two first terms are Cooper ones, whereas all two 
last terms are Peierls ones; so, we omit the spin depen- 
dence, which is given by the above equations, for each 
term. One gets 



dgo, . 



8A+4|c| v 



(E5o(c,±(A+isi) + i±£, T (A+4) + 4 £ )5o(^T(A+ i f i ) + ^ £ ,T(A+ i i i )-^ £ ) 



+ E 5to(c,±(A+^) + ^,T(A+^) + i±£)gt 7r (c, T (A+M) + i^, T (A+M)^i_£)") 
+ 8Aw(E9o(T(A+M )+ £±l !±(A+ M )+ £±i iP )g ( T( A + M )+ £^, T(A+ M)_£_l, p )) 
1E?r (E^(T(A+^^) + 4 i -A*/,±(A+i£±^i)+4i-Afc /) p) 

^(T(A + ^^) + ^ + A fc/ . T (A + ^M)-^+A fc/ , P+ 2A fe/ )) 



± 

A 

8A+4|p+2Afe/ 



8A$4W\ ( E ff/o(^±(A+^) + 4 £ ^(A+^) + i±£).g /0 (c,T(A+^) + ^£,T(A+M)„i 3 £) j 
+ 8A+4|c+2A M ( J 560( C ,±(A+^^) + i±£-A fe/ , T (cA+^M) + i± £ -A fc/ ) 

^( C+ 2A fc/ , T (A + ™^) + ^ + A fe/ , T (A+^±^)-^ + Afe/)) 
+ 8AT^H ( J 9f0(T( A +¥)+4 i >±(A+¥)+4 i .f)9/0(T(A+^)+^,T(A+^l)_^, p ) 
+ Effto(T(A+M ) + £^ :±(A+ M ) + £^, p )^( T(A+ J|i ) + £_l iT(A+ M ) _£ i d iP )") 



mo , , v 

( J 3o( c ' ± ( A+ ^ i ) + i i £ ^(A+l|i) + i±£)5 t o(c,T(A+^i) + ^£,T(A+lfi)-^£) 
+ E5to(c,±(A+iii) + ^,T(A+l|i) + ^)^(c, T (A+M ) + ^, T(A+ i£i ) _i_ £ )) 
+ 8A~^ ( E 5t0 (T ( A+ M ) + ^ ± ( A+ M ) + £^ 

+ E5/o(T(A+M )+ £+i, ± ( A+ M )+ £+i, p )5 to ( T(A+ M )+ £_i, T(A+ M ) _£ i i iP ) S ) 
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dgw 
dl 



(c,l,p) = 



8A+4|c| 



(E*o(c,±(A+l£i)+^+Afc / , T (A4^ 

^(ES6o(c,±(A + ^±^) + ^, T (A + ^^) + ^ 



8A+4|c+2Afc 



+ 8A4^(E5o(T(A+¥) + 4 A +Afc/,±(A+M) + 4i+Afc / , P ).g bo (T(A+M) + ^-Afc^ 
+ 8A+4^W(»°^ (A+ ^^ 



6. z couplings 



r 



Charge/Spin representation (a = C, S 1 ). Then, the spin 
dependence simply writes, for each one 



dz a 



dl 



9 a Zq 



Here are the RG equations for the couplings z, 
in (k,c,p) representation. The z sc couplings should 
be written in the singlet/triplet representation (a = 

s,t), and the z DW couplings should be written in the and will therefore be again omitted. One gets 



dzf c 



J"(C, k) = 3x ^ R (^Effo(c,±(A+lfi) + f-fe, T (A+^i) + f-fe)zS C ( c ,±(A + H) + S ) 
+ E.9to(c,±(A+M) + f _ fe ^ ( A+i|i) + f-fc)zS C (c,±(A+M) + S )^ 



sc 



-(c,fc) 



4A+2|c| 



Eff/o(c,±(A+ J | i ) + §-fe,T(A+lfi) + f-fc)4 C (c,±(A+J|i) + f) 



- 4A+2|c+2A fc/ | Egw(^(A+^)+|-A%-^(A+^) + f-At f -^ b °(c+2Afc / ,±(A + ^)+l + At J 



+ 



DW 



7^*0 = 



, 3x ^E3o(±(A+M) + f+fc , ± (A+^)-f-fc,p)z UW (p,±(A+^)-f) 



4A+2|p+2Afe/ 



E 56 7 r(±(A+^M) + f +A fe/+fe ,±(A+ 



|p + 2AJe f 



)-§+Afc / -fc,p+2Afc / )^ W (p+2Afc / ,±(A+ lP+ , r )-E-Ak f ) 



dz DW / 

-dT&V = 4Aw(E5/o(±(A+^) + f+^±(A+M)_ f _ fc , p ) z DW^ ±(A+ M ) _P) 
+ E5t7r(±(A+M) + | +fe , ± (A+i|i)-f-fe,p)^ W ( P ,±(A+M)_|)^ 



c. x couplings 



Here are the RG equations for the susceptibilities x, with the same spin dependence as the corresponding z couplings, 
which is again omitted, 



sc 



dl 



(«) = 4TO ( E 4 C ( 9 ,±(A+M ) + f )2 + ^ z SC (9 ±(A+ M ) + f )2) 
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dv sc 

-^(«) = -Tdm E 4 C («.±(A+¥) + i ) 2 - 4A +2 | g A + 2A fc/ | E Z S % + 2A fc/ , ±(A+ ^M) + f +A fc/ ) 2 



M W - _ 4A+2| 9 | E z {-iM^+W+t) - 4A+2|g-2Afc / | E Z tt (.-9+ 2 Afc/ ,±(A+ ^-) + f-Afc / j 



DW 



-(9) 



A 

4A+2|ij 



' ( E ^ W (-«,±(A+^) + f f + E Z DW (- g ,±(A + ^) + ^ 
V 4- + 



C. Symmetries 



a. Ordinary symmetries 



If we apply the conjugation symmetry C to the two- 
particle coupling S, we get : 



S(P[, Pi, P 2 , Pi) = S(-Pi, -P 2 , -P 2 , --Pi') , 



if we apply A, we get : 



Sc(P2,Pl,P 2 ',Pl') = 

-2Sc(Pi, P 2 , P 2 ', P{) - 3Ss(Pi, P 2 , P 2 ', P[) 

Ss(P 2 ,Pi,P 2 ',Pi') = 
Sc(Pi, P 2 , P 2 ', PD + 2Ss(Pi,P 2j P^,PD , 



if we apply ^4', we get : 



Sc(Pi,P 2 ,Pi,P^) = 
-2Sc(Pi, P2,P^,P{) - 3S s (Pi, P 2) P^, P[) 

S S (P 1 ,P 2 ,P 1 ',P^) = 
S c (Pi, P 2 , i^, PQ + 2gg(Pi 1 P2 1 jlPQ , 



and, finally, from parity P conservation, we get 

3(Pi,p 2 ,pLp{) = s(-Pi, -p 2 , -p^, -pO . 

Note that AA' simply gives 3(P 2 , Pi, P[, P^) = 

S(P u P 2 ,Pi,P{), 

For the SC instability coupling, we will write the two- 
dimensional interaction vector Q = (Qii,Q_l), and add 
a discrete variable 9 = 0, 7r, which indicates whether the 
R particle is on the 0-band (9 = 0) or the 7r-band (9 = 
7r); this way, one can distinguish 0-0, 0-7r, tt-0 or it-it 
processes (use Fig. [21 for help). 

If we apply P, we get 

Z sc (-Pi, -P 2) (-Q\\,Q±), 9) = Z sc (Pi, P 2 , (Q||, Qx), 0) , 

if we apply A or A 1 (note that in i/sc> the term with 
incoming momenta P\ and P 2 is conjugate to that with 
outgoing momenta P\ and P 2 ), we get 



^ C (P 2 ,P 1 ,(Q||,O),0) 
zf (s) (P 2 , P l5 (Q ||, f -Q ± ),0) 

zf (9) (p 2 , p l5 (g ,|,f -Q ± ),e) 
z t scfe) (p 2 , p l5 (g i,,f -Q ± ),fl) 

Z t SC(/H) (P 2 ,Pi,(Q||,f -Q_l),0) 



ZSC(P 1 ,P 2 ,(Q || ,o),e) a = s,t 

Z s s G ^\p u P 2 ,(Q h Q ± ), v -6) 

-Z S s C{9) (P 1 ,P 2 ,(Q 1] ,Qx),tt-9) 

Z t scfe) (P 1 ,P 2 ,(Q||,Q i _),7r-0) 

_ Z f^)(p 1 ,p 2 ,(Q || ,Q x ), 7r -^ 



Finally, it is interesting to note that Zf c (singlet) and 
Zf (triplet) change sign under S and are invariant under 
C, while Z( C and Z|p (both triplet) do the opposite. 

For the DW instability coupling, we use 
the same notation. Note that Q|| writes 

— (pi — P2 + %kfe) for intraband processes, and 

— (pi — p 2 + kfo + fc/7r) for intcrband ones. If we apply 
CS, we get : 

Z DW (-P 1; -P 2 , (2k fe - Q\\,Q±),9) = 
Z DW (P 1 ,P 2 ,(-2k fe + Q h Qx),9) , 



and if we apply AS, we get 

Z™(P 2 ,P 1 ,(2k f9 ~Q h Q 1 _),9) = 
±Z DW (P 1 ,P 2 ,(-2k fe + Q ll ,Qx),e) 

Z DW (P 2 , Pi, (k f0 + fc /7r - Q,|, f - Qx), 0) = 
TZ DW (Pi, P 2 , (-fc /0 - + Qy, Qx), 7T - 0) , 

where ± reads + for site ordering, and — for bond order- 
ing. 
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b. Supplementary symmetry 

When we apply the special symmetry C to two-particle 
couplings S, we get : 



S(fc/o + k f „ - P u kfo + k f7r - P 2 ,k f0 + k f „ - P' 2 , k f0 + k u - P[) = Q(P u P 2 ,Pl,, P[) 
When we apply the special symmetry C to SC instabilities Zsc, we get : 

Za C(r) (~fc/o -kf„- Pukfo + k fv -P 2 ,{-Q h Q),9) = ±Z s a c(r HP 1 ,P 2 , {Q\\,0),n-e) 

^a C(r) (-fc/0 - kf* - Pukfo + k f7T - P 2 , (-Q| h J - Q ± ),6) = ±Z S a G ^(P 1 ,P 2 , (Q|| , Q ± ), 7T - 9) 

where ± reads + for a = s, T = s or for a = t, T = p x , and — for a = s, T = d, g or for a = t, T = f x ,fy 
When we apply the special symmetry C to DW instabilities %DW: we get : 

Z DW (-fc /0 - k u - P u k f0 + k f „ - Pi, (2k fe - 0||, 0), 6) = ±Z DW (P U P 2 , (-2k f[v _ e] + Q h 0),tt - 0) 

Z DW (-k f0 - k fv - P u k f0 + k u - P 2 , {-kfo -k fv -Q h j- Q±_),9) = ±Z DW (P 1 ,P 2 , {-kf - k fn + Q h Q±),ir - 9) 
where ± reads + for site ordering, and — for bond ordering. 

c. gb orbits 

Here are the 8 first orbits of the coefficient, in (c,l,p) representation : 

{(0, 0, -2Ak f , 0), (0, -2Ak f , 0), (-2Afc/, 0, 0), (-2Afc/, -2Ak f , -2Afc/)} are sym. equiv. 

{(-4Afc/, -4Afc/, -2Ak f ), (-4Afc/, 2Ak f , 0), (2Afc/, -4Ak f , 0), (2Afc/, 2Afc/, -2Afc/)} id. 

{(0, -4Afc/, 2Afc/), (0, 2Afc /; -4Ak f ), {-2Ak f , -4Ak f , -4Ak f ), {-2Ak f , 2Ak f , 2Ak f )} id. 

{(-4Afc/, 0, 2Ah f ), (-4Afc /; -2Afc /7 -4Afc/), (2Afc/, 0, -4Afc/), {2Ak f , -2Ak f , 2Ak f )} id. 

{(0, -4Ak f , -2Ak f ), (0, 2Afc/,0), (-2Ak f , -4Ak f ,0), (-2Ak f , 2Ak f , -2Ak f )} id. 

{{-AAkf, 0, -2Ak f ), {-AAk } , -2Ak f , 0), (2Akf, 0, 0), (2Akf, -2Ak f , -2Ak f )} id. 

{(0, 0, 2Afc/), (0, -2Afe/, -4Afc/), {-2Ak f , 0, -4Afc / ), {-2Ak f , -2Ak f ,2Ak f )} id. 

{(-4Afc/, -4Afc /; 2Afc/), (-4AA;/, 2Ak f , -4Afc / ), (2Afc/, -4Afc /; -4Afc / ), (2Ak f , 2Ak f , 2Ak f )} id. 



D. Fourier Transform 

Creation and annihilation operators If one writes 
the creator of a particle of spin a, located in real space 
at position i (i e {1, • • • , N}), on chain j (j = ±1), the 
representation in the momentum space writes 

L p ,0,a = ^(-k fo +p,0)a = J2^ e ' iiP ~ kf0)ta 

ij 

^p,ir,a — M?(_fe /lr+Pi i ) CT — / JVija e 



J^p,TT,a — * (fe/7r+P,f)o' — / y J y»j<r e 

ij 

with the notations of the text. ^ k(T stands for the ab- 
solute momentum representation, while L and i? stand 
for the relative momentum representation. These rela- 
tions are given for annihilation operators, one must take 
the complex conjugation to obtain those for the creation 
operators. 

The reverse relations simply write, in terms of the ^ 
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operators, 

tv* = /" ^ eiMP (*(P,o)^+i*(P,f), CT ) ; 

but one can also express them in terms of the L and 7? 
operators. Then, one can check that this transformation 
is the inverse of the first one. 

Electron- electron pair operator In real space, the SC 
order parameters are the mean value of the electron- 
electron pair operator, which writes 

o a (X) = J2 ^ x >x' CT <r(x,x')T CT Q CT , . 

To each Q = (Q\\,Q±) corresponds a Fourier compo- 
nent 

O q (Q) = e~ QX Q (X) = e- i{aQ ^ i+hQxi/2) O a {X) . 

X ij 

We only keep the componants Q which lead to singu- 
larities; as explained in the main text, they are (0, 0) and 
(±Akf, So, using the short notation for the first 
and 7r± for the second, one gets 

o a (o) = J2° a ( ai > i) + ° Q K -i) 

i 

and 

O a (n±) = ^2 ~ie^ iAk f a (O a (ai, 1) - O a (ai, -1)) 

i 

where the main factor O a (ai, 1) ± O a (ai , — 1 ) = 
£w (±l)% iff ^ w r(a(i - i'),6(j - i')/2)0 is the 

mixed representation of the pair operator 45j . 

Eventually, the ipija can be expressed in terms of the 
^Pa, so that the componants write 

o a (Q)= e ^q(p)e*p.^q-p.-'^' 

a 0—0, 7T aa' 



where p = (P, |) and we will also use Q = (Q\\,Q±)- 
Be careful that, for instance, with Q = (Akf,j) and 
p = (p — fc/o,0), and thus ^ p . a = L p ,Q ai the calcula- 
tion of ^Q.p ^ is not immediate; one gets ^Q- P .a' = 

^ (Akf-p+k f0 ,f),a' — R2Ak f -p,-K,a' ■ 

With T(X, X') = Su'Sjji, one finds z (p) — 1 (singlet 
0-condensate of s symmetry), and z 7T± {p) = — i (sin- 
glet 7r-condensate of s symmetry). With T(X,X') = 
5ii'Sj,-j', one finds Zo(p) = cos(9) (singlet 0-condensate 
of d symmetry) and z w± (p) = i cos(#) (triplet 7T- 
condensate of f y symmetry). With T(X, X') = 
8i,i> T i5j,-ji , one finds z {p) = cos(aP) cos(9) (sin- 
glet 0-condensate of extended d symmetry), as well 
as z (p) = — isin(aP) cos(#) (triplet 0-condensate of 
f x symmetry), and z n± (p) = sin(a(P =p ^ X ))cos(6») 
(singlet 7r-condensate of g symmetry) or z„ ± (p) = 
ie ±I— J~ cos(o(P =F cos(9) (triplet 7r-condensate of 

extended f y symmetry). With 

T(X,X') = 8 i:i r T iSjji, one finds zo(p) = cos(aP) 
(singlet 0-condensate of extended s symmetry) or 
2o(p) = — isin(aP) (triplet 0-condensate of p x symme- 
try), and z v± (p) = — ie ±l— ^~ cos(a(P =F -^-)) (sin- 
glet 7r-condensate of extended s symmetry) or 2 Tr± (p) = 

; . Ak f a 

— e i~ sin(a(P =F — ^-)) (triplet 7r-condensate of p x 
symmetry) . 

Electron-hole pair operator It is almost the same, 
with the product of a creation and an annihilation op- 
erators, be careful, however, that, in reciprocal space, 
one gets : 



E L ^*p^Q +PCT '*(p)^< • 

aa' a 
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